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Abstract 

We investigate explicit functions that can produce truly random numbers. We 
use the analytical properties of the explicit functions to show that certain class 
of autonomous dynamical systems can generate random dynamics. This dynamics 
presents fundamental differences with the known chaotic systems. We present real 
physical systems that can produce this kind of random time-series. We report the 
results of real experiments with nonlinear circuits containing direct evidence for 
this new phenomenon. In particular, we show that a Josephson junction coupled 
to a chaotic circuit can generate unpredictable dynamics. Some applications are 
discussed. 
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1 Introduction 

Over the last three decades or so, a revolution has happened in the devel- 
opment of science. We are talking about Chaos theory [1-10]. In the chaotic 
regime, the behavior of a deterministic system appears random. This find- 
ing has forced many experimentalists to re-examine their data to determine 
whether some of the random behaviors attributed to noise are due to deter- 
ministic chaos instead. 
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Chaos theory has been successfully applied to many scientific and practical 
situations [1-10]. 

In the philosophical realm, however, the importance of this development was 
that chaos theory seemed to offer scientists the hope that almost "any" random 
behavior observed in Nature could be described using low-dimensional chaotic 
systems. Random-looking information gathered in the past (and shelved be- 
cause it was assumed to be too complicated) perhaps could now be explained 
in terms of simple laws. 

The known chaotic systems are not random [11]. If the previous values of a 
time-series determine the future values, then even if the dynamical behavior is 
chaotic, the future may, to some extent, be predicted from the behavior of past 
values that are similar to those of the present. The so-called "unpredictability" 
in the known chaotic systems is the result of the sensitive dependence on initial 
conditions. It is not an absolute unpredictability. 

Truly random systems are different from the chaotic ones. Past sequences of 
values of a random dynamical variable that are similar to present ones tell as 
much or little about the next value as about the next hundredth value. The 
so-called nonlinear forecasting methods for distinguishing chaos from random 
time-series are based on these ideas [11]. 

Recently, we have introduced explicit functions that produce truly random 
sequences [12-14]. For instance, let us define the function: 



where z is a real number and 6 is a parameter. 

For an integer z > 1, this is the solution to some chaotic maps [12-14] (see 
Fig. 1(a)). For a non-integer z, function (1) can produce truly unpredictable 
random sequences whose values are independent. 

Functions (1) with non-integer z cannot be expressed as a map of type 



In the present letter we address the following question: can an autonomous 
dynamical system with several variables produce a random dynamics similar 
to that of function (1)? We will present several dynamical systems with this 
kind of behavior. We will report the results of real experiments with nonlinear 
circuits, which contain direct evidence for this new phenomenon. We discuss 
some applications. 
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Fig. 1. First-return maps produced by function (1): (a) z = 5; (b) z = 7/3. 
2 Random functions 

Let us discuss first some properties of function (1). We will present here a 
short proof of the fact that the sequences generated by functions (1) are un- 
predictable from the previous values. This proof is presented here for the first 
time. However, a more detailed discussion of the properties of these functions 
(including statistical tests) can be found in Ref. [12-14]. 

Let z be a rational number expressed as z = p/q, where p and q are relative 
prime numbers. 

We are going to show that if we have m + 1 numbers generated by function 
(1): X , Xi,X 2 , X 3 , X m (m can be as large as we wish), then the next value 
X m+1 is still unpredictable. This is valid for any string of m + 1 numbers. 

Let us define the following family of sequences: 

X { n k ' m) = sin 2 

where k is an integer. The parameter k distinguishes the different sequences. 

For all sequences parametrized by k, the first m + 1 values are the same. This 
is so because 

X^ m) = sin 2 

for all n < m. Note that the number kp n q( m - n ) i s a n integer for n < m. So we 
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can have infinite sequences with the same first m+1 values. 



Nevertheless, the next value 
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is uncertain. 

In general, X^+i can take q different values. These q values can be as different 
as 0, 1/2, V2/2, 1/e, 1/tt, or 1. From the observation of the previous values 
X Q , Xi,X 2 , X 3 , ...,X m , there is no method for determining the next value. 

This result shows that for a given set of initial conditions, there exists always 
an infinite number of values of 9 that satisfy those initial conditions. The 
time-series produced for different values of 9 satisfying the initial conditions 
is different in most of the cases. Even if the initial conditions are exactly the 
same, the following values are completely different. This property is, in part, 
related to the fact that the equation sin 2 9 = a, where < a < 1, possesses 
infinite solutions for 9. 

We should stress that from the observation of a string of values X , X±, X 2 , X 3 , X, 
generated by function (1) it is impossible to determine which value of 9 was 
used. 

Figures 1(a) and 1(b) show the first-return maps for z = 5 and z = 7/3. 

For z irrational (we exclude the numbers of type z = m 1 ^), the numbers 
generated by function (1) are completely independent (see Fig. 2(a) that 
shows the first-return map for z — e). After any string of m + 1 numbers 
X , Xi, X 2 , X 3 , ...,X m , the next outcome X m+1 can take infinite different val- 
ues. 

The numbers produced by function (1) are random but are not distributed 
uniformly. The probability density behaves as P(X) ~ 1/ \j X{1 — X). If we 
need uniformly distributed random numbers, we should make the following 
transformation Y n — | arcsin \[X^. In this case P(Y) = const (see Fig. 2(b)). 

It is important to mention here that the argument of function (1) does not 
need to be exponential all the time, for n — > 00. In fact, a set of finite sequences 
(where each element-sequence is unpredictable, and the law for producing a 
new element-sequence cannot be obtained from the observations) can form an 
infinite unpredictable sequence. See the discussion in the following paragraph. 

So if we wish to produce random sequences of very long length, we can de- 
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Fig. 2. First-return maps produced by function (1): (a) z = e (irrational); (b) z = e, 



termine a new value of parameter 6 after a finite number N of values of X n . 
This procedure can be repeated the desired number of times. It is important 
to have a nonperiodic method for generating the new value of 9. For example, 
we can use the following method in order to change the parameter 9 after each 
set of N sequence values. Let us define 9 S = AW S) where W s is produced by 
a chaotic map of the form W s+ i = f(W s ); s is the order number of in a 
way that s = 1 corresponds to the 9 used for the first set of N values of X n , 
s = 2 for the second set, etc. The inequality A > 1 should hold to ensure the 
absolute unpredictability. In this case, from the observation of the values X n , 
it is impossible to determine the real value of 9. 

After a careful analysis of functions (1), we arrive at the preliminary conclusion 
that (to produce unpredictable dynamics) the main characteristics for any 
functions are the following: the function should be able to be re-written in the 
form 



where the argument function / (n) grows exponentially and the function h(y) 
should be finite and periodic. This result allows us to generalize this behavior 
to other functions as the following 



where P(t) is a periodic function. 

However, a more deep analysis shows that (to produce complex behavior) the 
function f(n) does not have to be exponential all the time, and function h(y) 




X n = h(f(n)), 



(6) 



X n = P(9z n ), 



(7) 
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does not have to be periodic. In fact, it is sufficient for function f(n) to be a 
finite nonperiodic oscillating function which possesses repeating intervals with 
finite exponential behavior. For instance, this can be a chaotic function. On the 
other hand, function h(y) should be noninvertible. In other words, it should 
have different maxima and minima in such a way that equation h(y) = a (for 
some specific interval of a, ot\ < a < a.2) possesses several solutions for y. 



3 Autonomous dynamical systems 



The following autonomous dynamical system can produce truly random dy- 
namics: 

J aX n , if X n < Q, , , 

Xn+1 -\bY n , ifX n >Q, (8) 

Y n+1 = cZ n , (9) 
Z n+l = sin 2 (nX n ) , (10) 



Here a > 1 can be an irrational number, b > 1, c > 1. We can note that for 
< X n < Q, the behavior of function Z n is exactly like that of function (1). 

For X n > Q the dynamics is reinjected to the region < X n < Q with a new 
initial condition. While X n is in the interval < X n < Q, the dynamics of Z n 
is unpredictable as it is function (1). Thus, the process of producing a new 
initial condition through the equation (9) is random. 

If the only observable is Z n , then it is impossible to predict the next values of 
this sequence using only the knowledge of the past values. 

An example of the dynamics produced by the dynamical system (8)-(10) is 
shown in Fig. 3. If we apply the nonlinear forecasting method analysis to a 
common chaotic system, then the prediction error increases with the number 
of time-steps into the future. On the other hand, when we apply this method 
to the time-series produced by system (8)-(10), the prediction error is inde- 
pendent of the time-steps into the future, as in the case of a random sequence. 
Other very strong methods[15], which allow to distinguish between chaos and 
random noise, produce the same result. 

Here we should make an important remark. Mathematical models can be of 
different types. For instance, natural phenomena can be described by differen- 
tial equations, difference equations, cellular automata, neural networks, etc. 
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Fig. 3. First-return maps produced by the dynamics of the dynamical system 
(8)-(10). (a) a = 7/3, b = 171, c = 1.5, Q = 1000/a; (b) a = e (irrational), 
b= 171, c = 1.5, Q = 1000/a. 

Even explicit functions can be mathematical models. If we consider func- 
tion (1) as the mathematical model of some dynamical process, then this 
process can be completely random. If all we know is the series of outcomes 
X , Xi,X 2 , then it does not matter how many values we already have, the 
next value is unpredictable. However, this is because from the observation 
of the values X , X 1: X 2 , X m , it is impossible to determine the "variable" 
Y n = 9nz n . This could be considered as a "hidden variable". Of course, we 
cannot say that this is just a problem of hidden variables, because not for any 
hidden variable Y n , the function X n = sin 2 (Y n ) is a random system. Here we 
have obtained necessary conditions for this phenomenon to occur. 

Let us extend this analysis to the dynamical system (8)-(10). In this case the 
completely random variable is Z n . The role of "hidden variable" is played by 
X n . If one could observe the series X n , then the complete randomness of the 
data set is lost. We could say that some projection of the complete set of 
variables onto a proper subset is necessary for the effect. 

It is important to notice that the dynamical system (8)-(10) is a well-posed set 
of difference equations with unique forward time evolution. Given the initial 
conditions (X , Y , Z ), the future evolution of the dynamics is fully deter- 
mined. However, if the only observable is Z n , then this variable will behave as 
a completely random time-series. How can one reconcile the unique forward 
time evolution with this randomness? The answer is related to the fact that 
only a subset of the variables are observed. 

With this result we are uncovering a new mechanism for generating random 
dynamics. This is a fundamental result because it is very important to under- 
stand different mechanisms by which the natural systems can produce truly 
random (not only chaotic) dynamics. 
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4 Pseudorandom number generators 

There is a large literature dedicated to pseudorandom number generators (see 
e.g. [16-38] and references quoted therein). A very fine theory has been devel- 
oped in this area and this theory has produced many important results. 

However, we should say here that the known pseudorandom number generators 
are not supposed to generate truly random numbers. 

In his important review article [17], James says : "Truly random numbers are 
unpredictable in advance and must be produced by a random physical process, 
such as radioactive decay" . 

In fact, pseudorandom numbers are produced using recurrence relations, and 
are therefore not truly random [17,18,24-26,33]. 

D'Souza et al [24] say in their paper: "Pseudorandom number generators are 
at best a practical substitute, and should be generally tested for the absence 
of undesired correlations" . 

Many known pseudorandom number generators are based on maps of type: 

X n +i = f (X n , X n _i..., X n _ r+ i) . (11) 

Now we will present the maps behind some of the most famous and best 
pseudorandom number generators. 

Multiplicative linear congruential generators [16,17] are defined by the follow- 
ing equation: 

X n +i = (aX n + c) modm (12) 

Some famous values for these parameters are the following: a = 23, m = 
10 8 + 1, c = 0; a = 65539, m = 2 29 , c = 0; a = 69069, m = 2 32 , c = 1; 
a = 16807, m = 2 31 - 1, c = 0; a = 1664525, m = 2 32 , c = (this is the best 
generator for m = 2 32 , according to the criteria of Knuth [20]). 

The Fibonacci-like generators obey the following equation: 

X n+1 = (X„_ p X„_ 9 ) modm, (13) 

where is some binary or logical operation. For instance, can be addition, 
subtraction or exclusive-or. 
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Other extended algorithms use equations as the following: 

X n+1 = (aX n + 6X„„i + c) mod m. (14) 

The add- and- carry generators are defined as: 

X n+ i = {X n - r ± X n _ s ± c) modm. (15) 

Among the high quality generators investigated in the famous paper [16] are 
the following 

X n+1 = (16807X n ) mod (2 31 - l) , (16) 
X n+ \ = (X n _ 103 .XOR.X n _25o) , (17) 
X n+ i = (X n _ W Q 3 .XOR.X n _i279) j (18) 

where .XOR. is the bitwise exclusive OR operator, 

X n = {X n _22 — A n _ 43 — c) , (19) 

here for X n > 0, c = 0, and for X n < 0, X n = X n + (2 32 - 5), c = I. 

All known generators (in some specific physical calculations) give rise to in- 
correct results because they deviate from randomness [16,24,25]. 

The problem is that these algorithms are predictable. 

An example of this can be found in the work of Ferrenberg et al [16]. They 
found that high quality pseudorandom number generators can yield incorrect 
answers due to subtle correlations between the generated numbers. 

Suppose we have an ideal generator for truly random numbers. In this case, no 
matter how many numbers we have generated, the value of the next number 
will be still unknown. That is, there is no way to write down a formula that 
will give the value of the next number in terms of the previous numbers, no 
matter how many numbers have been already generated. 

The authors of paper [16] related the errors in the simulations to the depen- 
dence in the generated numbers. Indeed, they are all based on maps of type 
(11). 

In the present paper we have shown that both the sequence of numbers X n 
defined by function (1) and the sequence of numbers Z n defined by the dy- 
namical system (8)-(10) cannot be expressed as a map of type (11). In fact, 
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Fig. 4. Noninvertible I-V characteristic. Two extrema. 

these numbers are unpredictable and the next value cannot be determined as 
a function of the previous values. 

Recently, simulations of different physical systems have become the "strongest" 
tests for pseudorandom number generators. Among these systems are the fol- 
lowing: the two-dimensional Ising model [16], ballistic deposition [24], and 
random walks [25]. 

Nogues et al [25] have found that using common pseudorandom number gen- 
erators, the produced random walks present symmetries, meaning that the 
generated numbers are not independent. 

On the other hand, the logarithmic plot of the mean distance versus the num- 
ber of steps N is not a straight line (as expected theoretically) after N > 10 5 
(in fact, it is a rapidly decaying function). 

D'Souza et al [24] use ballistic deposition to test the randomness of pseu- 
dorandom number generators. They found correlations in the pseudorandom 
numbers and strong coupling between the model and the generators (even 
generators that pass extensive statistical tests). 

One consequence of the Kardar-Parisi-Zhang theory is that the steady state 
behavior for the interface fluctuations (in ballistic deposition in one dimension) 
should resemble a random walk. Thus, a random walk again serves as a good 
test for pseudorandom numbers. 

We have produced random walks using the numbers generated by our systems. 
The produced random walks possess the correct properties, including the mean 
distance behavior (d 2 ) ~ N. 

The present paper is not about random number generators. In the present pa- 
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Fig. 5. Noninvertible I-V characteristic. Many extrema. 

per we discuss a new phenomenon: the fact that unperturbed physical systems 
can produce truly random dynamics. 

Of course, one of the applications of this phenomenon is random number 
generation. 

The art of random number generation requires more than the randomness of 
the generated numbers. It requires good programming skills and techniques 
to obtain the desired distributions for the numbers. 

The functions and systems described in this paper can be used to create very 
good random number generators. Algorithms designed for this purpose along 
with the statistical tests will be published elsewhere. 

In this section we only wished to present a theoretical comparison between 
the sequences produced by the pseudorandom-number-generator algorithms 
(11)-(19) and the systems described in the present paper. 



5 Experiments 

When the input is a normal chaotic signal and the system is an electronic 
circuit with the I-V characteristics shown in Figures 4 and 5, then the output 
will be a very complex signal. 

In Ref. [39] a theory of nonlinear circuits is presented. There we can find 
different methods to construct circuits with these I-V characteristic curves. 

The scheme of this composed system is shown in Fig. 6. A set of equations 
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Fig. 6. Scheme of a nonlinear system where a chaotic voltage source is used as the 
input signal for a nonlinear circuit with a noninvertible I-V characteristic. 

describing this dynamical system is the following: 

X n+1 = F l (X n ,Y n ), (20) 

Y n+1 = F 2 (X n ,Y n ), (21) 

Z n+1 = g(X n ), (22) 



where the equations (20) and (21) describe a normal chaotic dynamics where 
the variable X n presents intermittent intervals with a truncated exponential 
behavior and g (X n ) is a function with several maxima and minima as that 
shown in Fig. 5. 

Figures 7(a) and 7(b) show nonlinear circuits that can be used as the nonlinear 
system shown on the right of the scheme of Fig. 6. 

The system on the left of the scheme can be a chaotic circuit, e.g. the Chua's 
circuit [40]. 

We have constructed a circuit similar to the one shown in Fig. 7(a). We pro- 
duced chaotic time-series using a common nonlinear map and then we trans- 
formed them into analog signals using a converter. This analog signals were 
introduced as the voltage-input to the circuit shown in Fig. 7(a). Similar re- 
sults are obtained when we take the input signal from a chaotic electronic 
circuit. 

The set of equations that describes one of our experimental situations is the 
following: 

X n+1 = aX n [1 - 0(X n - q)} + bY n Q(X n - q), (23) 
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Fig. 7. Nonlinear circuits with noninvertible I-V characteristics, (a) The resistors 
possess R = 2.2fcf2, the source voltage in the battery is 1§V and the twin transistors 
are 2iV2222 with f3 = 140). The I-V characteristic of this circuit is shown in Fig. 4. 
(b) Another circuit with a similar I-V characteristic. 



where W n = ^ - 1, q = f , s = 10, b = 7, a = | , d = 3, 6(x) is the Heaviside 



The first-return maps of the sequence Z n produced by the theoretical model 
(23)-(25) and the experimental time-series produced by the nonlinear system 
of Figs. 6 and 7(a) are shown in Figures 8(a) and 8(b). 

When the nonlinear circuit has an I-V characteristic with many more maxima 
and minima, e.g. Fig. 5 (and this can be done in practice, see Ref. [39]), we 
can produce a much more complex dynamics. 

Nonlinear chaotic circuits can be described successfully by discrete maps as 
the equations (20)- (22) (see e.g. [41]). 

However, in some cases it can be very helpful to have a physical situation with 
a model based on a set of well-posed ordinary differential equations. 

In Ref. [42] we can find several models for chaotic circuits as the following: 




(24) 



z n+1 = ml - 3W r 



(25) 



function. 



dX 
~dt 



a[f(x)-Y], 



(26) 




(27) 
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Fig. 8. Modelling vs experiment, (a) Numerical simulation of the dynamical system 
(23)-(25); (b) first-return map produced with the real data (current measurements) 
from experiments using the scheme of Fig. 6, where the circuit of the right is the 
one of Fig. 7(a). 

f = /5r + 7 Z, (28) 
f =,(*). (29) 



N 

where / (x) = -X 3 + cX, g(X)= £ atf ■ 

i=i 

A comparison of different time-series is shown in Fig. 9. The fixed parameter 
values are: a = 285.714, (3 = 1499.25, c = 0.144, 7 = -0.51325. Note that 
when the I-V characteristic of the circuit shown on the right of Fig. 6 is a 
function with many extrema, the produced time-series is more complex (see 
Fig. 9(b)). 

Using our theoretical results we can make a very important prediction here. A 
nonlinear physical system constructed with chaotic circuits and a Josephson 
junction [43] can be an ideal experimental setup for the random dynamics that 
we are presenting here. 

It is well-known that the current in a Josephson junction may be written as 

I = I c sin 0, (30) 



where 



£ = *V ( 31 ) 
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Here is the phase and V is the voltage across the junction. Note that Na- 
ture has provided us with a phenomenon where the sine-function is intrinsic. 
Although we have already explained that other noninvertible functions can 
produce similar results, it is remarkable that we can use this very important 
physical system to investigate the real consequences of our results with func- 
tion (1). In a superconducting Josephson junction k is defined through the 
fundamental constants k = 2e/h. However, in the last decades there have 
been a wealth of experimental work dedicated to the creation of electronic 
analogs that can simulate the Josephson junction [45-48]. In that case k can 
be a parameter with different numerical values. 

We have performed real experiments with a nonlinear chaotic circuit coupled 
to an analog Josephson junction. 

In our experiments we have used the Josephson junction analog constructed by 
Magerlein [48] . This is a very accurate device that has been found very useful 
in many experiments for studying junction behavior in different circuits. The 
junction voltage is integrated using appropriate resetting circuitry to calculate 
the phase 0, and a current proportional to sin <p can be generated. The circuit 
diagram can be found in [48]. 

The parameter k is related to certain integrator time constant RC in the 
circuit. So we can change its value. This is important for our experiments. 
We need large values of k in order to increase the effective domain of the sine 
function. In other words, we need the argument of the sine function to take 
large values in a truncated exponential fashion. This allows us to have a very 
complex output signal. In our case the value of k is 10000. 

The voltage V(t) across the junction is not taken constant. This voltage will 
be produced by a chaotic system. In our case we selected the Chua's circuit 
[40] . For this, we have implemented the Chua's circuit following the recipe of 
Ref. [44]. 

The scheme of the Chua's circuit constructed by us can be found in Fig. 1 of 
Ref. [44]. 

The following components were used: C\ = 10 nF, C2 = 100 nF, L = 19 mH, 
and R is a 2.0 KQ trimpot. 

Chua's diode was built using a two-operational-amplifier configuration sug- 
gested in [44]. 

In our experiment, the voltage in C\ was used as the driving signal for the 
Josephson junction. We were interested in the famous double scroll attractor 
attained with R ss 1880 Q. 
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Fig. 9. Time-series generated by the dynamical sys- 
tem (26)-(29): (a) g (x) = 4X 3 - 3X; (b) 
9 ( x ) = 2k (88179X 11 - 230945X 9 + 218790X 7 - 90090X 5 + 15015X 3 - 693X). 

The differential equations that describe our experimental system are the fol- 
lowing: 

^=6.3(V 2 -Vi)-9/(Vi), (32) 

^T = 0.7(Vi-V 2 ) + J L , (33) 
at 

f = -7*. (34, 

£ = (35, 

^=sin0, (36) 



where /(K) = -0.5 [Vi + 0.3 (|V^ -h 1 1 — | Vi — 1 1 )] and fc = 10 4 . Notice that 
this system has been rewritten using adimensional variables (see Ref. [44]). 

The results of the experiments are shown in Fig. 10 which is the first-return 
map of the time-series data produced by direct measurements of the junction 
current. The time-intervals between measurements was 10 ms. This system can 
produce unpredictable dynamics. Figure 11 shows the results of the simulation 
of dynamical system (32)- (36). 
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Fig. 10. First-return map of the time-series generated with real data from an exper- 
iment with an analog Josephson junction coupled to the Chua's circuit. 



6 Conclusions 



In conclusion, we have shown that functions of type X n = P (9z n ), where P(t) 
is a periodic function and z is a noninteger number, can produce completely 
random numbers. Certain class of autonomous dynamical systems can gener- 
ate a similar dynamics. This dynamics presents fundamental differences with 
the known chaotic systems. We have presented real nonlinear systems that 
can produce this kind of random time-series. We have reported the results of 
real experiments with nonlinear circuits containing direct evidence supporting 
this phenomenon. 

Besides the fundamental importance of these findings, these systems possess 
many practical applications. For instance, game theory tells us that in certain 
competitive situations the optimal strategy is a random behavior. Specifically, 
it is necessary to limit the competition's ability to predict our decisions. We 
can produce randomness using the discussed systems. Another example is 
secure communications [49]. In this context, the most important application 
of our systems is masking messages using random signals [50]. In some cases, 
when we use the usual chaotic systems, the messages can be cracked because 
the time-series are not truly unpredictable. 

Now we will analyze very general ideas. 

Just to facilitate our discussion (because it is always important to have a 
name), we will call the phenomenon studied in this paper "deterministic ran- 
domness" . The words "deterministic randomness" have been used (metaphor- 
ically) in the past as a name for chaos. However, the known chaotic systems 
are not random. So we think this is a good name for the present phenomenon. 
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Fig. 11. First-return map of the time-series generated by the dynamical system 
(32)-(36). Here /= 

Deterministic randomness imposes fundamental limits on prediction, but it 
also suggests that there could exist causal relationships where none were pre- 
viously suspected. 

Deterministic randomness demonstrates that a system can have the most com- 
plicated behavior that emerges as a consequence of simple, nonlinear interac- 
tion of only a few effective degrees of freedom. 

On one hand, deterministic randomness implies that if there is a phenomenon 
in the World (whose mechanism from first principles is not known) described 
by a dynamical system of type (8)-(10) or (20)-(22), and the only observable 
is a physical variable as Z n , then the law of this phenomenon cannot be learnt 
from the experimental data, or the observations. And, situations in which the 
fundamental law should be inferred from the observations alone have not been 
uncommon in physics. 

On the other hand, the fact that this random dynamics is produced by a rel- 
atively simple, well-defined autonomous dynamical system implies that many 
random phenomena could be more predictable than have been thought. 

Suppose there is a system thought to be completely random. From the ob- 
servation of some single variable, scientists cannot obtain the generation law. 
However, suppose that in some cases, studying the deep connections of the 
phenomenon, we can deduce a dynamical system of type (8)-(10) or (20)-(22). 
In these cases, some prediction is possible. 

In any case, what is certain at this point is that some dynamical systems can 
generate randomness on their own without the need for any external random 
input. 
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